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1 Introduction 



A lot of theoretical and experimental effort has recently been made towards 
understanding the adsorption on substrates equipped with geometric and 
chemical structure. This research exposed a variety of phenomena and asso- 
ciated problems, to mention the filling transitions [1-4], position and droplet's 
size-dependent contact angles and interfacial morphology structures [5-10], 
discontinuous changes in droplet shapes as function of volume [11-12], or the 
infiuence of substrate's structure on the existence and order of wetting tran- 
sitions as compared to the homogeneous and planar substrate case [13-17]. 
In many such systems the interfacial morphology leads to the existence of 
three-phase contact lines and in consequence to the line tension [TH]. This 
quantity's properties, and especially its behaviour close to the wetting tran- 
sitions stimulated many discussions and controversies in recent years [19-29]. 

In this paper we are concerned with semi-infinite fluid in a thermody- 
namic state close to its bulk liquid-vapour coexistence, in the presence of a 
substrate consisting of a single stripe-like inhomogeneity placed on an other- 
wise chemically homogeneous and planar solid surface, see Fig.l. The system 
is translationally invariant along the stripe, say in the y-direction, and the 
equilibrium adsorption morphology varies only in the x-direction perpendic- 
ular to the stripe, the width of which is denoted by 2R. Thus the chemical 
inhomogeneity of the substrate imposes nonuniformity of the liquid-like layer 
adsorbed at the substrate. Its thickness varies in the x-direction and induces 
linear contribution to the fluid free energy. This contribution evaluated per 
unit length of the stripe will be denoted 77, and called the line tension. It 
is a function of temperature T, the inhomogeneity width 2i?, and addition- 
ally depends on quantities characterizing both the stripe and the rest of the 
substrate. These additional quantities will be represented by the relevant 
wetting temperatures of two parts of the substrate. The y = const section 
of the system can be looked upon as a two-dimensional droplet sessile on 
a one-dimensional substrate inhomogeneity, see Fig. 2. We are particularly 
interested in contributions to the droplet line tension resulting from its flnite 
size - their size and temperature dependence. 

In addition to analyzing the droplet morphology and the contributions to 
the line tension we perform geometric constructions of the relevant contact 
angles. These constructions correspond to different choices of the three phase 
contact line and lead to different predictions concerning the behaviour of the 
contact angle in the limit i? — > 00. In effect, only one of them turns out to 
have practical meaning. 

This work is arranged as follows. In Section 2 we specify the model and 



discuss its limitations. A set of constrained equilibrium profiles parametrized 
by the droplet's height at its centre is constructed and the constrained excess 
line free energy expressions as functions of droplets' height, inhomogeneity 
width 2i?, temperature, and the corresponding wetting temperatures are 
derived. In Section 3 we construct the expansion for equilibrium droplet 
height in the limit of large R. We recover some of the earlier specific results 
|8fl| together with new expressions for the corrections up to the order which 
is necessary to calculate the leading i?-dependent contributions to the line 
tension. In Section 4 the expressions for the line tension are derived and 
discussed. In Section 5 we analyze the contact angles obtained from the 
equilibrium liquid - vapour interfacial shapes. 

2 The model 



We consider a planar, chemically inhomogeneous substrate remaining in con- 
tact with a fluid. The system's thermodynamic state is inflnitisimally close 
to the fluid's bulk liquid- vapour coexistence, i.e., [i = fij^t^ where fisat denotes 
the chemical potential at liquid-vapour coexistence. The inhomogeneity has 
the form of a stripe of width 2R. The substrate is translationally invariant 
in the ^/-direction defined by the inhomogeneity's orientation, see Fig.l. The 
wetting temperature of the stripe T^i is lower than that of the remaining 
part of the substrate Tw2, i-e., Twi < Tw2, so that the liquid phase is pref- 
erentially adsorbed on the stripe. The system's temperature T is assumed 
lower than Tw2 so that the equilibrium thickness of the adsorbed layer l{x) 
converges in the limit x -^ ±oo to a finite value determined by the properties 
of substrate 2 alone, see Fig. 2. 

A suitable description of this system, valid on length scales much larger 
than the bulk correlation length ^b is provided by the interfacial Hamiltonian 
model ISH, 132] 



H[l] = dx dy 



-{vir+v{i,x) 



(2.1) 



Here l{x, y) denotes the interface position, a - the liquid-vapour surface ten- 
sion, and V{l^x) - the effective interfacial potential. The borderline between 
substrates of types 1 and 2 is assumed to be of the order of ^b and the 
crossover in V{l^x) is anticipated to occur over a distance of the same mag- 
nitude for systems with short-ranged intermolecular forces |j2i}i, iMj- Conse- 
quently, we model V{l,x) as 

V{1, x) = Q{R - \x\)Vi{l) + e(|x| - R)V2{1), (2.2) 



Type 1 substrate 




Type 2 substrate 



Figure 1: A schematic illustration of a planar substrate containing a stripe- 
like inhomogeneity. The wetting temperature of the stripe (T^i) is assumed 
lower than the substrate 2 wetting temperature {Tw2)- The domain width 
(2i?) is considered much larger than bulk correlation lengths. 




Figure 2: Schematic plot of the y = const section of the fluid density in 
the vicinity of the stripe, and the corresponding equilibrium shape of the 
liquid- vapor interfacial proflle denoted by l{x). 



where Vi(/) and V2(/) are the effective interfacial potentials corresponding 
to homogeneous substrates of type 1 and 2, respectively, and Q{x) is the 
Heaviside step function. For short-range forces each of the potentials Vi{l) 
{i = 1,2) has the exponential form (see e.g., [11], |30] ) 



Vi{l) = r,e~'/«^ + 6e-2'/«s + 



(2.3) 



where Ti = a{T — T^^j), a being a positive constant, and ^b denotes the bulk 
correlation length in the adsorbed liquid phase. The positive parameter b 
controls repulsion of the interface from the substrate at short distances. For 
simplicity, within our model, it has the same value for substrates of both 
kinds. 

Thorough this work we apply the mean-field approach according to which the 
equilibrium grand cannonical average value of the adsorbed layer's thickness 
< l{x) > is identified with l{x) which minimizes the Hamiltonian (j2.1|) . One 
should also note that the full drumhead expression in the interfacial Hamil- 
tonian (see e.g., |^) is substituted with the gradient term in p. HI . This 



approximation is valid for small differences between both substrates which 
in turn is controlled by the difference ri — T2. The equilibrium profile T{x, y) 
is invariant with respect to translation in the y direction and in addition 
l{x,y) = l{x) = l{—x). From now on the considered values of x are limited 
to X G [0, 00 [. 

Minimizing the functional 1)2.1^ leads to the Euler-Lagrange equation 

dH _ dV 
dx"^ dl 
together with the derivative continuity condition 



(2.4) 



dl 
dx 



=R- 



dl 
dx 



(2.5) 



--R+ 



The minimization procedure is performed in two steps. First we solve Eq. ()2.4|l 
subject to the boundary conditions: ^\x=o = 0, lima-^oo ^(a;) = Itt2, where 
It,2 corresponds to the minimum of V^il). The droplet's height at the centre 
Iq = l{x = Q) is for the time being kept as fixed but arbitrary parameter. This 
way we obtain a set of constrained equilibrium profiles denoted by {/(x, /q)}; 
they are parametrized by /q. In the second step the equilibrium droplet's 
height Iq and thus l{x) = l{x, Jq) are determined by minimizing H[l{x, Iq)] 
with respect to Iq. 



For this purpose two cases must be distinguished: 1) Vi(/o) < 0; and 
2) Vi{lo) > 0. In case 1), upon applying ()2.2|l . p.3p we obtain the following 
formulae for the constrained equilibrium profiles corresponding to fixed value 
lo- 



l{x,lo) 




2yi 

2b 

T'2 



^(n + (ri + 26e"^"/«s) cosh( 



-2Vi(lo) 



?£ 



■ik _i 



^2 



{^-c)/iL 



whilst for case 2) we have 



is 




1 



2Vi{Iq) 

T2 



ri + (ri + 26e-'«/^s)cos( 



2Vifo) 



T2 



ix-C)/!iB 



))] 


for X < i? 


for X > i? 


(2.6) 


for X < i? 




for X > R. 

(2.7) 



The above expressions are valid for x > and negative x-values are reached 
via /(— X, /o) = /(x, /o). The constant C is in each case determined by the 
continuity condition of /(x, Iq) at x = R. Note, that /(x = R~; Iq) is uniquely 
given by Iq. We require that /q is such that 



l{x = R-,lo) >i 



7r2 



(2.8) 



The necessity to make this additional though rather natural assumption is a 
consequence of the discontinuous form of the interfacial potential V{l,x), 
see Eq. ()2.2jl . which allowed us to solve Eq. ()2.4ll ior x < R and x > R 
independently. The profiles l{x, Iq) are decreasing functions of x, concave 
for X < R and convex for x > R. 
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Figure 3: An example of constrained equilibrium droplet profile l{x,lo). The 
plot parametres were chosen as follows: 2b/a = 1, Ti/a = 0.1, T2/a = —0.01. 
The lengths x, Iq = 7.2, R = 100 are expressed in units of ^b- The profile's 
first derivative is discontinuous at x = R (not visible in this figure's scale). 
The magnitude of the derivative discontinuity decreases to zero for Iq —■ Iq. 



At this stage one may calculate the constrained line contribution to the 
free energy per unit length in the ^/-direction as function of /q and R; it is 
defined as 



H{lo,R) 



dx 



a fdl{x, Iq) 
dx 



+V{xJ)-Vi{Ui)e{R-x)-V2{K2)Q{x-R) 

(2.9) 

Substituting expression ()2.6|] for case 1) and ()2.7jl for case 2) we obtain the 
following formulae 



H{k,R) 



V2^ 



^2 I 2Vi(/o) \ , .R 

2b n + in + 26e-'o/?fl) cosh(a)y ''i^'o;^^^ 



+ 2t-iW— arthf 



2b 



b , ic , a 

-V^i(/o) ^2 



-2aVi(/o) 



sinha 



(2.10) 




r2log 



rfR 



» (n + (n + 26e-W..) eosh(a))j + ^i^e(-n) 
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for case 1), and 






+ H^bM^Wo)'"'''''^) "^ V^^^)^r^- + 




-^2 



46^1 (/o) 



ri + [n 



26e~'»/^^ 



sma 

+2fc;-'o/eB +cos« 

(2.11) 

46 ^B 



for case 2), where a = \/\ — ^-^|^. Note that the set of physically relevant 
values of Iq is limited by the condition p.Sp . which means that Iq > lomin, 
where /omm is given by l{x = R,lQmin) = 1-k2- The inhomogeneity width R 
is arbitrary, nevertheless in further calculations we shall assume ^||i <^ i?, 
where ^||j denotes the mean-field interfacial correlation length corresponding 
to substrate of type i = 1,2, ^,7^ = V-'{lT,i)/a. 

The formulae ^^, ^J^, (l2Tnil . and (irTlTl form the starting point for the 
subsequent analysis aiming at the determination of the droplets' height, line 
tension, and the contact angles. The line tension rj as function of Tj and R is 
given as the minimum of if (/q, R) with respect to Iq. The corresponding value 
of drolet's height is denoted by Iq. Consequently, the equilibrium interfacial 
shape l{x) = l{x, Iq), and r] = H(Iq, R). 



3 The equilibrium droplet height 



The scaling behaviour of the equilibrium droplet height Iq for i? — > cxd and 
Ti > 0, was obtained by Parry et al using conformal properties of equation 
()2.4p in the case 6 = [30]. However, this technique fails when the term 
~ e~'^^I^B jjj Vx{p) must be taken into account. As we show below, the term 
neglected in |2Qj contributes to the corrections to the leading behaviour of Iq 
which are essential to determine the character of the asymptotic expansion 
of line tension r] for large R. In the following subsections we expand the 
functions — g^°' ' around the properly chosen values /q of /q for different 

temperature ranges. We assume ^o = /q + 5/o, (^^o ^ ^o ? ^^^ solve the 
appropriate equation for bl^, thus determining Iq up to terms necessary for 
the calculation of the leading correction to r] due to finite values of R. 
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3.1 The Ti > regime 



In this case the numerical analysis of the constrained line tension H{lo, R) 
given by Eq. (l2.11|) shows that the equilibrium droplet height fulfils the con- 
dition a{lo{R)) — i> 7r^, as i? -^ oo. The range of physically relevant values 
of Iq is limited by the condition ()2.8p . It follows that only Iq > lomin should 
be considered, where /omm is defined by l{x = R,lomin) = ^2- We note that 
H{Iq, R) diverges at Iq 



^0 < hmin and one can check that asymptotically 



/* 



Omin 



and simultanously /q 



1*0, for R 



oo. One may not perform the 



dH{lo,R) 
dlo 

this value is smaller than Iq for finite R (see Fig. 4). 



expansion of "'^g^''"^ around the value of Iq which corresponds to a 



Tc, as 



H(l^R) 



k) lomin 1 



lo 



Figure 4: A sketch of H{Iq, R) for fixed R. The physically meaningful values 
of droplet height fulfil Iq > lomin, where /omm is determined by the condition 
l{x = RJomin) = K2- The function H{Iq,R) diverges at Iq < lomin- 



Instead, one expands — i,°' ' around /, 



dlo 



(0) 




/g, which fulfils 



m, 



n + (ri + 2be-^o'/^B^ cos[a{li°\ R)] = 0. 
The solution of Eq. (|3.1jl takes - for ^^ ^ 1 - the following form 



-CVcb 



an^es 



2(tK| 



^V2ab^B I, 2/oV 



Oi- 



(3.1) 



(3.2) 



2rii?2 V Ti R ' 2^" " ' ' rf R"^ ' ' WiR' 

dlo 

up to terms of the order {^f^Y- In this way one 
obtains the equilibrium droplet height 

2c2 , ,^ ,^ _t 



Similarly, by expanding ""^0'"-> around the above value l'^' , we solve the 
equation — g°' ' 



-lo/iL 



2rii?2 ^, ^ R^ R2 



0{{ 



TiR' 



(3.3) 



where ki = ^^^806^^^^^, ^^2 = 8abC 



B 



8-r 



+ 



- 2Tf +5t^ T2-2T1 T^ +Ti 



T1T2 Tir|(2Tf-3nT2+T|) 



If in Eq. (|3.3jl one neglects the correction terms of the order 
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;?!)' ^^d 



~ (^^)^, then one recovers the result obtained by Parry et al via conformal 
transformation |3fl| . 



The droplet height /q = ^b 



log 



1^2 



2TlRf_\ Kl K.2-2l^l 



TT^a^l I R R? 



exhibits log- 



arythmic divergence for _R — > oo. The first correction (~ 1/-R) is positive. 
Both Ki and ^2 depend on h and thus cannot be obtained within the confor- 
mal scheme [201 • Note, that the asymptotic expansion was performed while 
keeping Tj {i = 1, 2) constant. Taking the limits ^b/R -^ and Xj/cr — > are 
non-commutative operations and - consequently - the above expansion (j^l^^jl 
is valid only provided ^^ <^ 1. This remark remains also relevant for the 
cases Ti = and ti < discussed below. 

3.2 The Ti = regime 

The plot of the function H{lo, R) for ti = case is qualitatively similar 
to the Ti > case, see Fig. 4. However, in the present case the quantity 
Iq corresponds to ailo^R)) — ^ |~ for i? — >^ oo. Following, we may expand 

H{lo, R) around Iq which fulfils the condition a(/g , R) = tc/2. After solving 
the equation — g°' ' = we obtain 

where 71 = — 2^||2, 72 = 4^m2- The explicit expression for Iq has the form 

9(t 2f ^ 

^^^ + R — W + •••]• The leading correction term, alike 
case Ti > 0, is positive. 

3.3 The Ti < regime 

Our first observation concerning the case ti < is that in the limit i? — > oo, 
the equilibrium droplet height lo{R) converges to a finite value /,ri correspond- 
ing to the minimum of the potential Vi(/). Indeed, as the inhomogeneity size 
becomes infinite, the height of the drop at its centre is fully determined by 
the properties of substrate type 1. The quantity Iq = It,i is therefore a 
natural choice of the value around which — q°' ' can be expanded for large 
R. Note that unlike the former cases ti > 0, the effective potential Vi{Iq) 
becomes negative close to Iq = lo{R) — l-jri and so the formula in Eq. ()2.in| l 
should be used instead of (12. mi . Solving — q°' ' = up to the relevant order 
leads to 

k = Li + 5ie-^/«iii + 52e-2«/«Hi + O(e-3^/«ii0 , (3.5) 



where 6i = 2^^^, and 62 = ^Sf. The leading correction is negative and 
equal for the special case ti = T2 as is expected on physical ground. 

To sum up the results obtained for the morphology of the liquid drop 
we note that the form of the To expansion for large R depends crucially on 
the value of the reduced temperature ti. For ri > which corresponds 
to complete wetting of type 1 substrate one has corrections proportional to 
powers of 1/R while for ri < we obtain corrections proportional to powers 



4 The line tension 



The expressions for the line tension ri{R) = Hijo, R) are obtained by substi- 
tuting the formulae for Iq into equations (I2.1()jl , and (J2.1H1 . The calculations 
are rather cumbersome and we refrain from qouting the expressions for rj in 
terms of the parametres of the expansion of lo{R) around R = 00. We only 
note that the leading /^-independent term in rj may be expressed in terms 
of the parametres ki, 71 or 61 for all cases ri > 0, ri = 0, ri < 0, respec- 
tively. The next-to leading corrections in /q modify r] at the level of leading 
i?-dependent terms. 

We obtain the following expressions: 



v>{R,n,r2) 




Tilog 



n -T2 



< 



B 



T"2log 



2(r2 - n) 



T-2 



''CB^2^ 



yi+°'<f» 



(4.1) 



for Ti > 0, 



r]o{R, Ti, T2) = -^bTs J — log(2) - a^B 



r^Hi 



8 i? V r' 



(4.2) 



for Ti = 0, 
and finally 



r]<{R,Ti,T2) =^B 



CbTi 




^1^0g( ^o- ^ ) +T-2l0g(^ 



n + T2 



2b I Ti + Ta 



2ri / ^ "V 2t2 



+ 



(4.3) 
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for ri < 0. 

The dominant, i?-independent terms in the above formulae may be ob- 
tained independently by considering an inhomogeneous substrate consisting 
of two homogeneous half-planes characterized by ri, r2, meeting at a; = 0. 
After calculating the corresponding excess linear free energy one obtains the 
dominant, i?-independent term. It is a non-negative function of Ti and T2, 
which vanishes for Ti = T2, see Fig. 5. 




-0.5 



0.5 



Ti/(T 



Figure 5: Plot of the leading i?-independent contribution to line tension 
denoted by r]^^^ as function of ri at fixed T2. The plot parametres were 
chosen as follows: a/2b = —T2/2b = 1. Note that for Ti = T2 one has r]^^^ = 
as expected on physical ground. The range of ri is restricted to ri > T2 in 
accordance with the analysis in the text. 



The central observation concerning the above formulae is the change in 
the character of the corrections to the "free" line tension 77*^°^ as the temper- 
ature Twi is crossed. The interfacial correlation length ^||i = ^r^ sets the 
length scale controlling the next-to-leading terms in ?7<, i.e., for ri < 0. Its 
divergence at ri = is accompanied by the appearance of the algebraic terms 
in the expansion of rj about 1/R = for ri > . The finite R corrections 
for Ti > are universal (at least up to ~ 1/R order) in the sense that they 
do not depend on the substrates' characteristics and are fully determined by 
a. Thus the divergence of ^||i at ri = is accompanied by the change of the 
expansion parameter of 77 from e'^/^H^ to ^b/R- 

The corrections are positive for ri < and negative for Ti > 0. They may 
be interpreted as the asymptotic term of an effective interaction potential 
between two regions of inhomogeneity of /(x), corresponding to x ~ i? and 
X R. 
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5 Remarks on the contact angles 



At the macroscopic level the contact angle of a liquid drop placed on a pla- 
nar solid surface is related to the interfacial tensions between the relevant 
phases. Accordingly, it can be evaluated from the Young's equation [TH]. At 
the mesoscopic level the 3-phase contact line is not defined unambigously and 
consequently the geometric construction of the contact angle is not unique. 
Moreover, generalizations of Young's equation to the case of chemically in- 
homogeneous substrates assume implicitely that the typical size of boundary 
regions between homogeneous domains is much larger than bulk correlation 
length C,B (see Ref.0 for detailed discussion), which is opposite to the case 
considered in this paper. 

There are at least two approaches towards the geometric definition of the 
contact angle for the system considered in this paper. The first definition 
exploits the interfacial profile's characteristics in the region of inhomogeneity 
(i.e., for |x| ~ R) while the second probes the profile's features at the centre 
of the droplet (i.e., x ~ 0). 

Within the first approach the equilibrium contact angle is defined via the 
profile's slope at the three phase contact line situated at a choosen position 
Xk- A natural - though not unique - choice for the considered system is 
Xk = R- Following, the contact angle 6 is obtained from 



dlix) 
tan ^ = ^ ' 



dx 



(5.1) 



x=R 



In the second approach one fits the circle y{x) = \J'R? — x^ + y^, tangent 
to l{x) at X = 0, and of the same curvature at x = as that of l{x)\ 
2^' l^^o = —^ = ^^2 |a;=o. The position of the contact line is then defined 
as the line of intersection of y{x) with the profile's asymptote at It,2, see 
Fig. 6. 

Within the present analytical approach the above two constructions can 
be achieved only in the case of large inhomogeneity width -R/^||i ^ 1. For 
the first choice one obtains 

for Ti < 

for Ti > ' ^^-^^ 

where 9i and 02 are contact angles for homogeneous substrates of type 1 and 
2, respectively, and for small, negative Ti one has 9i ~ |rj|/v2a6, i = 1,2. 
It is obvious on physical ground that the second construction leads for -R/^||i ^ 
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Figure 6: Phenomenological construction of the contact line and the contact 
angle. A circle y{x) fitted to the equilibrium droplet profile at a; = intersects 
the asymptote /^2 at x'^ defining the three phase contact line position. The 
contact angle 9' is defined as the corresponding angle of intersection. 



1 to the vanishing contact angle. This is confirmed by the analytical expres- 
sions 



6' 



2|rik/(/. 



T1+T2 2ab 
^ v/log[-(2/(a&))l/2r2fi/^] 
V2 R 



U^z^^e-i^/% forn<0 



n 



y^°g[^ 



fl2. 



R 



for Ti = 



for Ti > 



(5.3) 



which display change of the convergence type to the value 9'{R = oo) = 

from exponential for ri < to ^^'^-^ — for ri > 0. In other words, the dis- 
crepancy between our predictions for 9 and 9' may be attributed to different 
behaviour of the corresponding positions Xk and x'^ of the contact line for 
R ^ oo. In particular, the difference Xk — x'j^ is of the order e^/*^^^iii^ for 
Ti < and ^/\ogR for ri > 0. 



6 Summary 



This article is concerned with morphology and line tension of equilibrium 
droplets placed on a planar substrate with a chemical, stripe-shaped inho- 
mogeneity. 

Applying the effective Hamiltonian approach, we first constructed a set 
of constrained equilibrium states parametrized by the droplet's height. The 
global minimum was then extracted from this set using perturbative expan- 
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sion valid in the limit of large droplet radii. The line tension expressions were 
evaluated and the character of the size-dependent correction shown to exhibit 
crossover from exponential to algebraic decay for increasing inhomogeneity's 
width, occuring at the strip's wetting temperature. This corresponds to di- 
vergence of the interfacial correlation length characterizing the wetted part 
of the surface, which acts as the length scale controlling line tension for low 
temperatures. 

In addition, on the basis of the obtained equilibrium liquid-vapour profiles 
l{x) we calculated the droplet contact angle on the heterogeneous substrate 
defined via tangents of l{x) evaluated at x = R. 

The main conclusions are restricted to the case of short-ranged forces and 
refer to the asymptotic regime -R/^||i -^ oo, and the thermodynamic state 
close to the fiuid two-phase coexistence (/i = figat)- ^^^ predictions are as 
follows: 

• The line tension rj related to the substrate's inhomogeneity is a function 
of the stripe's width R. For large R the character of the leading in- 
dependent term depends on the temperature. At low temperatures, 
i.e., for T < T^i the interfacial correlation length ^||i represents the 
characteristic length scale controlling the i?-dependent contributions to 
line tension, which are of the type e'^'^n^. For T > Twi the correlation 
length ^||i is infinite and the i?-dependent contributions to 77 become 
algebraic and universal, i.e. independent of substrates' properties, and 
fully determined by the liquid-vapour interfacial tension a. 

• The i?-dependent contributions to r] may be interpreted as the effec- 
tive interaction potential between the interfacial profile inhomogeneity 
regions corresponding to droplet boundaries at a; ~ i? and x ~ —R. 
At T = Twi the leading i?-dependent correction to 77 changes sign 
from positive in the low temperature regime T < Ty^i to negative for 
T > T]vi. This means that the interaction represented by these terms is 
repulsive as long as T does not exceed T^i and attractive for T > T^i. 
The interaction magnitude is governed by the inverse of the correlation 
length 1/C||i = |ti|/v2ct6 for T < Twi and by the interfacial tension a 
in the case of T > Twi- 

• In order to characterize the interfacial morphology in terms of the con- 
tact angle we analysed two definitions of the contact angle. They turn 
out to be highly sensitive to the choice of the three phase contact line. 
A particular choice of this line positioned at the domain boundary leads 
to a contact angle expressible in terms of contact angles on homoge- 
neous substrates. Another choice, inspired by experimental procedures 
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and based on droplets' properties at their centres, leads - in the present 
context - to contact angles asymptotically (for R -^ oo) equal for all 
temperatures. 
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